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摘 要 : 计算 机 对 人 类 情绪 与 情感 的 识别 研究 已 经 成 为 了 脑 机 接口 领域 的 研究 热点 。 通 过 分 析 人 类 在 生活 中 的 各 种 情 


感 状 态 ， 提 取 脑 电信 号 的 特征 并 对 情感 状态 进行 识别 、 分 类 是 情感 智能 化 领域 的 重要 方向 。 针 对 基于 音乐 视频 诱导 的 
情感 数据 集 DEAP 进行 了 研究 ， 提 取 脑 电信 号 的 频 域 特征 后 ， 提 出 了 采用 加 速 近邻 梯度 算法 《APG) 和 正 交 匹配 算法 
(OMP) 求解 稀疏 编码 的 稀疏 表示 分 类 模型 进行 情感 分 类 ,并 与 支持 向 量 机 算法 SVM) 做 效果 比较 。 实 验 结果 表明 ， 


APG 算法 通过 Ll 范 数 正则 近似 求解 以 其 快速 的 收 你 速度 在 情感 数据 集 上 有 着 较 好 的 分 类 表现 ， 而 OMP 算法 与 SVM 
算法 的 分 类 效果 相差 无 几 ， 实 现 了 情感 脑 电 信号 的 分 类 。 
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Classifying emotional EEG using sparse representation method 
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Abstract: Computer recognition of human emotion has become a hot topic in the field of brain computer interface (BCD in 


recently years. By analyzing the various emotional states in people’s life, extracting the features of EEG and classifying 
emotional states is an important direction in the field of emotional intelligence. Based on the emotion data set induced by the 
music video, this research extracted the frequency-domain features of EEG. After that, the Accelerated Proximal Gradient (APG) 
and Orthogonal Matching Pursuit (OMP) algorithms for the sparse representation method were adopted to classify the EEG 
signals. By comparing with other algorithms, the experimental results show that the APG with L1 norm performs well in the 
emotion data set with fast convergence speed, and the greedy idea based OMP algorithm can achieve the same effect with other 


algorithms. The comparative analysis in this paper show the effectiveness and feasibility of the proposed method for emotional 


EEG signals classification. 
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知 科学 以 及 人 机 自然 交互 关键 的 一 环 。 根 据 现 代 认 知 科学 与 心 
理学 ， 情 感 是 人 脑 对 外 界 的 客观 事物 的 一 种 反映 ， 它 是 人 对 客 
在 过 去 很 长 的 一 段 时 间 内 ， 人 类 情感 的 研究 一 直 被 认 知 科 ” 观 事物 是 否 满足 自己 的 需求 而 产生 的 态度 体验 。 因 此 ， 研 究 人 
学 研究 者 所 忽略 。 直到 20 世纪 末 , 情感 作为 认 知 过 程 中 重要 的 赂 与 情感 的 关系 对 于 未 来 计算 机 识别 人 类 情感 的 研究 至 关 重 要 。 
一 部 分 才 被 人 们 所 重视 。 与 人 类 运动 的 研究 相 比 ， 研 究 参与 情 章 乐 是 人 类 思想 感情 和 艺术 创作 的 结合 ， 从 十 到 今 无 数 的 
感 活动 的 神经 元 的 具体 功能 难度 更 大 。 神 经 生理 学 和 心理 学 的 。 音乐 家 用 音乐 诠释 各 式 各 样 的 情感 ， 可 想 而 知音 乐 对 人 类 情感 
研究 发 展 在 一 定 程度 上 把 情感 的 由 来 从 神秘 莫 测 中 揭露 出 来 。 的 产生 与 诱发 是 有 着 特殊 意义 的 。 脑 电波 信号 (EEG) 在 很 大 程 
情感 是 人 类 特有 的 生理 活动 ， 包 括 情感 的 表达 、 情 感 的 识别 以 ” 度 上 能 够 反映 人 的 心理 和 生理 状态 的 改变 。 通 过 研究 人 脑 中 枢 
及 情感 的 转变 等 。 计 算 机 对 于 人 类 情感 的 识别 是 人 工 智 能 、 认 ”神经 系统 产生 的 生物 电信 号 发 现 ， 人 在 情绪 波动 或 者 受到 不 同 
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形式 的 感官 刺激 时 ， 脑 电波 信号 会 产生 特定 模式 的 变化 。 脑 电 
波 信号 是 一 种 自发 的 、 非 平稳 、 非 线性 的 神经 电 活动 ， 其 频率 
变动 处 于 低频 节律 与 超 低 频 节律 之 间 ， 可 划分 为 多 个 波段 ， 常 
见 的 波段 频率 有 : delta 波 (0.5~3 Hz)、theta 波 (4~7 Hz)、alpha 波 
(8~13 Hz)、beta 波 (14~30 Hz)，gamma 波 (31~47 Hz) 等 。 幅 值 范 
用 为 5~300 4 A。Karthick 等 人 上 通过 采集 聆听 古典 音乐 和 摇滚 
音乐 以 及 安静 闭 眼 休息 三 种 状态 下 的 额 叶 脑 电 信号 ， 采 用 去 趋 
势 波动 分 析 (detrended fluctuation analysis, DFA) 算 法 和 多 尺度 焙 
(multiscale entropy MSE) 分 析 脑 电信 号 发 现 有 、 无 音乐 时 脑 电 
言 号 有 着 显著 的 差异 .Nie 等 人 外 用 不 同情 感 类 别 的 奥斯卡 经 典 
电影 镜头 片段 ， 设 计 了 一 个 以 电影 片段 作为 刺激 的 实验 ， 通 过 
提取 被 试 脑 电信 号 频 域 特征 ， 开 进行 正 负 情绪 的 分 类 。Sander 
等 人 BI 以 40 个 一 分 钟 时 长 的 音乐 视频 (music videos, MV) 作 为 
刺激 素材 ， 采 集 脑 电信 号 和 生理 信号 进行 情感 识别 研究 ， 并 且 
记录 了 被 试 面 对 不 同 的 音乐 视频 刺激 时 的 主观 评价 以 及 被 试 在 
实验 时 的 面部 表情 视频 ， 建 立 了 情感 分 析 数 据 库 “DEAP”。 
目前 脑 电信 号 的 采集 设备 主要 是 脑 电 图 仪 ， 它 可 以 记录 放 
在 头皮 的 电极 所 检测 出 来 的 微弱 脑 电 信号 ， 通 过 电极 导 联 耦 
合 到 差 动 放大 器 进行 适当 的 放大 并 且 数 字 化 ， 最 后 通过 与 其 配 
套 的 PC 上 的 记录 系统 记录 下 信号 数据 ,关于 头皮 电极 的 位 置 ， 
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DEAP 数据 集 


图 2 效 价 -唤醒 度 二 维 空间 


DEAP (database for emotion analysis using physiological 
signals ) 数据 库 是 Sander KoelstraB3] 等 人 利用 多 模 态 的 生理 信号 
和 音乐 视频 进行 情感 分 析 的 数据 集 。 其 中 生理 信号 包括 被 试 的 


多 通 


感 状态 分 类 方法 。 本 文采 用 


道 脑 电信 号 CEEG) 和 外 围 生 理 信 号 。 这 是 一 个 公开 的 情 
感 类 EEG 数据 集 , 研究 人 员 可 以 基于 此 数据 库 来 测试 各 种 的 情 


了 此 数据 库 中 的 部 分 数据 进行 研究 。 


实验 选取 了 32 名 被 试 ， 男 女 各 一 半 ， 平 均 年 龄 26.9 岁 。 


在 实验 前 告知 被 试 实验 的 


的 以 及 实验 的 各 个 细节 。 实 验 在 有 


可 控 照 明 的 实验 室内 进行 ， 


试 在 实验 时 情绪 不 受 打扰 .被 试 坐 在 距离 屏幕 将 近 1 米 的 位 置 ， 


有 许多 放置 法 如 Montreal 法 ，Cohn 法 以 及 Gibbs 法 等 。 但 应 用 
最 多 的 是 10/20 系统 法 ， 即 国际 脑 电 图 学 会 建议 采用 的 标准 电 
极 安放 法 。 采 集 脑 电 信号 的 电极 数 通常 有 16 导 、32 导 、64 导 、 
128 导 。 在 图 1 中 展示 了 64 导电 极 的 分 布 情况 。 


Re 


对 


图 1 64 导电 极 安放 示意 


1 ”情绪 的 维度 


情感 往往 产生 于 特定 的 环境 当中 ， 而 人 的 情感 又 是 变化 无 
常 的 。 目 前 研究 中 最 常用 是 基于 情感 的 认 知 理论 的 二 维 空间 模 
型 : 效 价 (valence〉 和 唤醒 度 〈arousal) 。 其 中 效 价 是 指 人 的 
情感 的 愉快 程度 ， 从 消极 情绪 到 积极 情绪 ;而 唤醒 度 指 的 是 人 
的 兴奋 程度 , 从 无 聊 到 兴奋 .丹麦 心理 学 家 卡尔 格 (Carl Lange) 
提出 喜欢 -兴奋 二 维 空间 内 ， 与 效 价 与 唤醒 都 二 维 空 间 模型 相似 ， 
横 坐 标 表 示 愉 快 程度 ， 纵 坐标 表示 兴奋 程度 ， 不 同 的 情绪 分 布 
在 二 维 空间 的 不 同位 置 ， 如 图 2 所 示 。 


实验 中 的 音乐 视频 在 17 
过 与 被 试 沟通 将 音乐 音量 
的 投入 到 实验 环境 当中 。 准 
开始 播放 40 首 音乐 视频 素材 。 


实验 室内 保持 安静 的 状态 以 保证 被 


英寸 的 液晶 显示 屏 中 播放 。 实 验 前 通 
周 试 到 恰当 大 小 ， 确 保 被 试 能 够 更 好 


E 备 工作 完善 之 后 ， 通 过 下 面 的 流程 


实验 流程 如 下 : a 显 示 2s 的 音乐 编号 ,以 提醒 被 试 当前 实 
验 进度 ; b)5 s 的 基线 时 间 ; c) 播 放 1 min 的 音乐 视频 ， 并 且 记 


录 生 


里 信号 ， 采 样 频率 为 512 Hz; d) 被 试 自我 评估 。 


其 中 被 试 进行 自我 评估 是 根据 自我 评估 系统 SAM (self- 
assessment manikins ) 表格 进行 的 ， 该 系统 记录 了 用 户 对 40 首 
音乐 视频 的 唤醒 度 ， 效 价 、 喜 欢 程度 和 控制 力度 的 评分 ， 每 项 


是 32 个 导 联 。 实验 数 ] 
率 降 至 128 Hz。 为 了 去 
滤波 至 4~45 Hz， 然 后 


根据 国际 10-20 


评分 的 评分 标准 均 为 1~9 分 。 


极 法 


， 实 验 用 于 采集 实验 信号 的 电极 帆 


居 已 经 经 过 了 简单 的 预 处 理 , 其 中 采样 频 
除 噪 声 ， 首 先 用 带 通 滤波 器 将 脑 电信 号 


] 盲 源 分 离 法 去 除了 眼 电 伪 迹 。 每 首 音 


乐 视频 对 应 的 脑 电 信号 被 截取 为 63 s， 其 中 60 s 为 音乐 视频 的 


播放 时 间 ，3s 为 基线 数据 时 间 。63 s 的 数据 共计 8064 个 采样 
点 ， 本文 在 处 理 数据 时 ， 去 除了 前 3 s 的 基线 数据 ,因此 采样 点 


降 为 7680 个 。 


鉴于 不 同 的 被 试 对 于 音乐 的 敏感 程度 不 一 样 ， 本 文通 过 分 


析 被 试 的 自我 评估 系统 表 (SAM) ,选取 了 10 个 效 价 和 控制 力 


度 评估 较 高 的 被 试 作为 测试 样本 , 他 们 分 别 是 : S04, S06, S09， 


S17, S19, S20,， S21, 


额 叶 皮层 在 情感 调节 和 意识 经 验 中 起 着 主要 作用 ， 故 本 文选 取 


了 处 于 前 额 的 FP1, FP2 电极 与 额 区 F3，F4 电极 的 四 个 导 联 做 


S22,，S26，S27。 


Davinson 和 Petrantonakis 等 人 6 9 的 研究 表明 ， 人 脑 的 前 


据 形式 


综 上 所 示 , 在 本 文 的 研究 当中 , 每 一 个 被 试 的 EEG 信号 数 


为 40x4x7680， 其 中 40 代表 有 40 首 音乐 视频 ，4 代表 


采样 信号 。 


original signal 
T 


10 位 被 试 的 数据 均 是 在 
3 所 示 为 S06 被 试 在 FP1 通道 


所 选 的 电极 ， 分 别 是 FP1，FP2，,F3，F4 电极 ; 7680 代表 被 试 
在 1 min 音乐 视频 下 的 脑 电 采样 信号 。 
MATLAB2014 下 处 理 的 ， 如 图 
的 原始 采 


上 1 h 
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a 5000 
图 3 S06 被 试 在 FP1 通道 下 的 采样 信号 
3 ”算法 流程 
3.1 算法 框架 
本 文 的 算法 流程 框架 如 图 4 所 示 。 


DEAP 情感 数据 集 


beta、gamma 五 1 
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3.2 ”z-score 数据 标准 化 


数据 标 ; 
来 说 ， 往 外 


同 的 ， 


z-score 数据 祭 准 化 


利用 小 波 变换 提取 delta、theta、alpha、 
波段 节律 


利用 快速 传 里 叶 变 换 计 算 各 个 节律 的 功 


光环 这 :| 县 
率 谱 密度 


分 类 对 比 


4 ”算法 流程 


舍 化 是 数据 处 理 前 的 


因 


项 基本 工作 ， 


据 过 SVM 算法 ， 稀 下 表 示 分 类 中 的 OMP 算法 、 
APG 算法 进行 


对 于 脑 电信 号 


时 分 析 结 果 的 影响 ， 


具有 多 个 评价 指标 ， 但 是 多 个 评价 指标 的 量 纲 是 不 
此 ， 为 了 消除 量 纲 对 于 数 提 


需要 对 


方法 


进 


本 文选 ] 


处 理 前 
标准 化 


行 标准 化 的 处 理 。 基 于 数 


统计 学 ， 


min-max 标准 


取 标 准 分 数 (z-score ) 对 数 ] 


化 、1log 函数 转换 、z-score 
昌 进 行 处 理 。 如 


放 上 
中 让 


的 数据 标准 
标准 
式 (1) 所 示 ， 


化 
化 [7] 等 。 


的 数据 为 x，/ 为 变量 的 均 


EE 为 村 


后 的 数据 。 被 处 理 后 的 脑 


忆 信 号 采 档 


FF 点 数 直 


k 准 差 值 ，z 值 


表示 
昌 被 缩放 到 [- 


1,1] 内 


3.3 


EEG 


时 频 分 析 方 法 , 它 
年 提出 的 。1987 各 


o 


小 波 变换 
言 号 是 典型 的 非 平稳 时 变 、 非 高 斯 信号 ， 而 传统 的 频 
域 分 析 方 法 要 求 信 号 是 平稳 信号 ， 因 此 ， 研 究 人 员 想 到 采 / 
频 分 析 方 法 。 小 波 变 换 (wavelettransform, WT) 是 一 种 典型 的 
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欣 ， 等 : 基于 稀疏 表示 的 脑 电 《〈EEG) 情感 分 类 
本 (1) 


时 


的 数 


E 论 引入 小 波 分 市 


FE， 法 国 


学 形式 是 有 Grossman 和 Morletltsl 在 1984 
到 信号 处 理 专 家 Mallet 将 多 尺度 分 析 的 


Ff 中 ， 其 中 包括 了 小 波 函 数 的 分 解 以 及 信 


波 函 数 重 构 等 ， 完 善 了 小 波 分 析 的 理论 体系 。 


情感 脑 


电信 号 是 大 脑 中 各 种 凶 


经 元 之 间 相 互 作用 的 随机 组 


组 合 的 非 线性 导致 了 脑 


电信 号 具有 随机 性 、 非 线性 、 非 平 


a 
口 


导 处 理 方法 为 傅 里 时 变换 ， 它 需要 利用 


/区 \， 


具有 一 定 的 局 限 性 。 小 波 变换 在 传 里 叶 


小 
和 频率 窗 都 可 变 的 分 析 方 法 ,更 


了 时 间 旬 


号 变化 较 慢 的 
氏 频 的 情感 脑 电 


些 非 平稳 信号 。 对 于 活动 频 
波 变换 可 以 通 
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电 信号 


表现 脑 电 的 在 小 尺度 下 


2 
人 
KU 
a 


口 


表现 高 频 


同时 呈现 脑 电 信 入 系数 和 近 


这 样 就 可 以 在 不 同 尺度 


上 提取 脑 电 不 同 节 律 的 波 


里 情感 脑 


小 波 变 


响 (WT) 


下 六 
I 


号 的 过 程 中 ， 小 波 变换 体现 了 对 信号 的 


的 定义 为 式 (2) ， 其 中 a 是 伸缩 系数 ， 


z 是 平移 系数 。 小 波 变换 对 信和 号 分 析 处 理 时 通过 母 小 波 函数 h(t) 


来 完成 的 ， 
变换 可 以 选 


其 中 Ce ) 表示 
择 出 被 分 析 的 信号 , 信号 通过 a 尺度 的 伸缩 .平移 ， 


其 共 斩 复 数 ， 母 小 波 函 数 在 时 域 中 的 


对 于 频率 而 言 也 是 相似 的 。 


1 im tt 1 
WT.m0= | "On d= Val xadn 过 0 


小 波 变换 通过 在 特定 范围 子 集 上 对 信号 缩放 和 平移 变换 ， 
将 原始 情感 脑 电信 号 分 解 到 不 同 尺 度 的 近似 分 量 和 细节 分 量 上 ， 
每 一 层 的 近似 系数 再 次 分 解 为 新 的 近似 分 量 和 细节 分 量 ， 以 此 
类 推 ,将 信号 分 解 成 n 层 ,每 层 对 应 一 组 近似 分 量 和 细节 分 量 。 


如 式 (3) 所 示 ， 对 于 


其 中 : 


大 个 细 


128HZ, 根据 奈 奎 # 
即 64Hz。 根 ] 
节 分 量 和 
Hz)、alpha 
Hz) 五 个 波 


原始 信号 x(t), 小 波 分 解 过 程 为 


N= C02 "D+ DD p02 6) 


Cg 为 第 n 


= A = 


节 分 量 ，o(t) 表 示 伸 缩 函 
本 文中 


层 的 第 个 近 


PF DEAP 情感 


以 分 量 ，Dj, 代表 第 j 层 的 第 


数 ，9(D 是 伸缩 函数 。 
尚 电 数 据 已 经 过 预 处 理 ， 采 样 频率 为 


居 式 (3) 


特 采 样 定理 ， 


近似 分 量 中 提取 出 : 
波 (8~13 Hz)、beta 波 (14~30 Hz)，gamma 波 (31~47 
段 ， 本文 选取 db4 作为 小 波 基 , 图 


有 效 频率 为 采样 频率 的 一 半 ， 
言 号 分 解 为 四 层 ， 从 每 层 的 细 
delta 波 (0.5~3 Hz)、theta 波 (4~7 


， 将 脑 电 


5 为 其 对 应 的 4 层 
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小 波 频 域 分 解 树 。 图 6、7 为 被 试 S06 原始 脑 电信 号 小 波 分 解 ” 随 着 时 间 移 动 的 时 域 信号 转变 为 功率 随时 间 变 化 的 频 域 信号 的 
对 应 的 近似 系数 和 细节 系数 。 谱 图 。 经 典 的 频谱 分 析 方 法 是 以 有 限时 间 段 数据 的 傅 里 叶 变 换 
为 基础 的 周期 法 ， 同 时 信号 采用 窗口 截断 对 信号 功率 谱 进 行 估 

计 。 根 据 维 纳 - 辛 钦定 理 ， 信 和 号 x(n) 的 功率 谱 和 它 的 自 相关 函数 


ER 可 构成 一 对 传 里 叶 变换 ， 它 的 定义 为 式 (4) ， 该 方法 通过 将 随 
机 信号 x(n) 的 N 个 观测 数据 看 做 是 能 量 有 序 的 序列 ， 计 算 x(n) 
的 傅 里 叶 变 换 , 然 后 取 其 幅 值 的 平方 除 以 N, 作 为 序列 x(n) 的 


SRC2(0~8HZ) SRD2(8~16HZ) 实 功 率 谱 估计 。 
SRC1(0~4HZ) SRD1(4~8HZ) I,00) 坟 |X CW) 2 (4) 


图 1 4 层 小 波 分 解 树 


x(o) 的 傅 里 叶 变换 具有 周期 性 ， 所 以 RN(w) 也 有 周期 性 ， 并 
如 图 5 所 示 ， 原 始 的 64HZ 信号 被 分 解 成 四 层 ， 每 层 包括 。“ 且 是 有 偏 估计 。 先 计算 N 个 数据 点 的 傅 里 叶 变换 〈 即 频谱 ) ， 
一 组 近似 分 量 SRC 与 细节 分 量 SRD, 其 中 alpha 波 分 布 于 细节 ”如 式 (5) 所 示 。 


分 量 SRD2 上 ，beta 波段 分 布 在 细节 分 量 SRD3 上 ，theta 波 分 AI ee 

、 X= xWe (5) 
布 于 细节 分 量 SRD1 上 ，gamma 波 分 布 于 细节 分 量 SRD4 上 ， n=0 
delta 波 分 布 于 近似 分 量 SRC1 上 。 然后 频谱 域 与 其 共 恩 函数 的 相 乘 , 得 到 其 功率 谱 , 如 式 (6) 


所 示 。 


Detail Coefficients 


| 2 
| 1 2 1 | i 
Ope Heeb wth DL 已 .oo 一 六 IO = *Ve (0) 
100 200 300 400 500 0 
吕 基于 经 典 的 谱 估计 方法 ， 本 文 分 别 计算 出 已 提取 出 的 五 种 
9 
8 ED J 
节律 波形 的 功率 谱 ， 如 图 ， 并 将 其 选取 为 特征 ， 四 个 通道 总 共 
5 T T T 
要 4x5=20 个 特征 。 
SA 
0 500 1000 1500 2000 
T T T I The PSD of delta 
已 0 上 
-1 时 上 二 二 025 10 20 30 40 50 60 70 
0 500 1000 1500 2000 2500 3000 3500 4000 The PSD of theta 
sampled signal i ] 
0 
网 i 
6 细节 系数 Es 2 hebsp of Sie 洲 多 
10 | 
oo 10 20 30 40 50 60 70 
The PSD of beta 
Approximation coefficient 4 r 
5 T T T T 2 | 
025 10 20 30 40 50 60 70 
可 0 | The SD of mide 
: 0 ds 
-0 100 200 300 400 500 "0 10 20 30 40 50 60 70 
10 + T T T ferquency/HZ 
Be 图 8 五 个 频带 的 功率 谱 
-10 1 1 1 1 
0 200 400 600 800 1000 -一 ， 
10 3.5 ” 稀 玻 表示 (sparse representation ) 
operehhwhedome thoests fete AAA 
0 元 i 二 7 本 原理 是 将 输入 的 信号 表示 为 从 字典 中 选取 的 过 完备 基 上 的 线 
5 T T T T T T T > 一 、 
性 组 合 ,其 中 对 于 一 个 N 维 向 量 x, 如 果 其 中 大 部 分 元 素 均 为 0， 
,0 
只 有 少 部 分 元 素 是 非 零 的 , 那 x 被 称 之 为 是 稀疏 的 。 Donoho 于 
5 1 | | 1 | 六 _ 的 
0 500 1000 i 3000 3500 4000 2006 年 提 出 了 压缩 感知 (CS) 理 论 使 言 号 的 稀疏 表示 得 以 发 


展 00。 压 缩 感知 使 得 信号 的 采集 和 测量 方法 不 再 受 奈 码 斯 特 
CNyqyist) 采样 定理 的 限制 ， 对 信号 的 带宽 不 再 严格 要 求 ， 取 


图 7 近似 系数 


3.4 功率 谱 分 析 而 代 之 的 是 数据 的 稀疏 性 5。 
情感 脑 电 信号 有 着 其 幅 值 小 、 噪 声 干 扰 大 的 特点 ， 因 此 ， 不 妨 把 一 个 数据 集 D 看 做 是 一 个 矩阵 , 若 将 其 每 一 行 看 做 


如 何 从 中 提取 有 用 的 特征 是 一 个 很 重要 的 问题 。 目 前 基于 频 域 ”一 个 样本 ， 每 一 列 看 做 是 一 个 特征 ， 那 么 特征 选择 过 程 中 要 考 
的 分 析 是 研究 脑 电 信号 和 临床 应 用 的 主要 分 析 方 法 。 通 过 小 波 。 虑 的 就 是 这 些 特征 的 “ 稀 疏 性 ”， 若 去 掉 矩 阵 中 与 学 习 任务 无 
变换 已 提取 出 情感 脑 电 信号 的 频 域 信息 。 频 域 分 析 可 以 将 幅 值 ” 关 的 特征 , 则 训练 过 程 的 难度 以 及 计算 和 存储 的 开销 也 会 减 小 。 
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另 一 方面 , 若 D 中 存在 很 多 零 元 素 , 但 它们 不 是 以 整 行 或 者 整 
列 的 形式 出 现 ,这 时 对 于 数据 集 来 说 , 它 同 样 是 具有 稀 疏 性 的 ， 
当 样本 具有 这 样 的 稀 玻 表达 形式 时 ， 不 仅 可 以 使 得 大 多 数 的 问 


b) 规 则 化 矩阵 4 的 列 ， 


使 其 具有 


0 解决 Ll 范式 最 小 化 问题 : 


minllxl|h subject to 


题 变 得 线性 可 分 ， 而 且 存 储 的 负担 大 大 降低 ， 因 为 稀 疏 矩阵 已 


量 ,表示 重 构 误差 限制 ; 


有 很 多 高 效 的 存储 方法 。 基 于 压缩 感知 理论 的 信号 采样 只 需 使 
用 不 到 原 信号 一 半 的 样本 数 ， 即 可 高 质量 地 重建 原 信 号 。 本 研 
完 利 用 情感 脑 电信 号 不 同 节律 的 稀 玻 特性 《或 可 压缩 性 ) 的 先 
验 条 件 ， 通 过 不 同 的 解码 算法 可 以 很 高 的 概率 重建 原始 脑 电信 
号 或 者 高 维 信号 ， 对 于 重要 的 节律 特征 波形 能 够 很 好 的 保持 
且 找 与 情感 高 度 相关 的 节律 特征 。 

信号 的 稀 玻 表示 问题 可 以 通过 求 
现 ， 如 式 (7) 所 示 。 


解 稀 玻 正则 化 问题 得 以 实 


Mings(®)E f(a)+ AdQ(a) (7) 


其 中 : 8 为 稀疏 正则 项 , 4 为 非 负 a 为 4 中 的 
向 量 。 
系数 a 具备 稀疏 性 ，Q 可 以 取 Lo 范 数 ， 如 式 〈8) 所 示 。 
Qa)= el 全 |{xlax 尖 0}| (8) 
其 中 |,| 为 集合 的 势 。 
f(a) 表 示 信 号 的 误差 项 , i 


参数 ，4 为 凸 集 ， 


通常 取 其 平方 误差 , 式 (9) 所 示 。 


f(a) = 了 |* 一 Do (9) 


稀 玻 表示 在 用 于 分 类 时 的 模型 被 称 之 
(Sparse Representation based Classification, SRC) [3 。 
主要 是 将 输入 信 
式 (10) 所 示 。 


为 是 稀疏 分 类 方法 
该 模型 
号 表示 为 稀 朴 的 过 完备 基 原 子 的 线性 组 合 ， 如 


min|| x llo subject to Ax=y 


其 中 : xj 表示 x 中 非 零 的 元 素 个 数 , 即 Zo 范 数 。 通 过 Zo 范 


数 表示 稀 疏 性 虽然 理论 可 行 ， 但 是 在 数学 上 很 难 实现 。 常 见 的 
做 法 是 将 Lo 范 数 转 换 为 L1 规 范 化 处 理 ， 如 式 (11) 所 示 。 


min||xll subjectto Ax=y (11) 

可 见 ， 式 11) 转换 为 一 个 最 优 凸 近似 问题 ， 和 型 的 凸 松 
弛 测量 是 通过 声 正 则 近似 得 到 ， 如 式 〈12) 所 示 。 

min g(a) 51x pal +4lalh (12) 

式 (12) 是 统计 学 中 的 Lasso 问题 ， 可 以 看 做 是 一 个 正则 

最 小 二 乘 问题 。 本 文采 用 原始 情感 脑 电信 号 作为 过 完备 基 ， 利 

用 过 完备 基 的 宛 余 特性 可 以 更 好 的 捕 提 信 号 脑 电 信号 中 的 本 质 


特征 这 一 特点 , 提出 了 采用 SRC 分 类 算法 对 情感 脑 电 信号 进行 
分 类 。 
SRC 分 类 框架 的 具体 流程 如 下 : 


a) 输 入 训练 样本 矩阵 4=[4 ,4 4]e 形 ”"” ,测试 样本 


yeR” ; 


dd) 计算 残 差 : 


6) 根据 重 构 样 本 和 


1(y)=||y—AS Wl, 12 


Ax 一 yl a, 其 中 为 党 


测试 样本 的 见 余 误差 确定 测试 样 


签 : identity(y) = arg min,(x(y)) 


稀 琉 正则 化 求解 通常 


了 两 种 方法 ， 


一 种 是 通 


解 近似 解 ， 另 一 种 是 将 其 转换 为 凸 优化 问题 。 


3.5.1 贪 焚 算 法 求解 


其 高 效 的 特性 在 压缩 感知 和 稀 玻 信号 重建 中 得 
多 次 迭代 字典 中 与 


人 
到 ) 2 应 上 43 巴 多 和 本 原理 是 通过 
相似 的 原子 来 允 近 原始 信号 。 贪 


禁 算 法 主要 包括 匹配 追踪 


(matching pursuit, MP )03 
pursuit，OMP ) 69 ， 


matching pursuit) 05] 等 。 


、 正 交 匹 配 追 踪 (Corthogonal matching 
分 级 正 交 匹配 追踪 〈stagewise orthogonal 


MP 算法 在 每 次 迭代 中 ,选择 与 当前 误差 向 量 
求 出 残 差 ， 然 后 选择 与 残 差 最 匹配 的 原 
L 体 步骤 如 下 : 


最 相关 的 向 量 里 ， 
影 到 选择 的 原子 上 ， 构 成 该 原子 的 


a) 初 始 化 误差 向 量 


r=x— Da 


系数 增 量 。 


b) 搜 索 与 误差 向 量 相关 性 最 大 的 原子 索引 : 


1 :一 arg max |d; x 


is[[1: | 


co) 更 新 重 构 系数 与 误差 向 量 : 
ay :一 ay 十 


gf1:=1 十 1] 若 t1<k ， 返 回 


1 
Lt:=% ;一 (d; 


. 


Fd 汉 


dr fl 
| di ; 
步 又 b)， 


否则 退 


子 投 


1 二 XX ， 连 代 次 数 1:==1; 


出 。 


MP 算法 中 ， 


其 他 原子 正 交 ， 


rd 二 0， 但 是 却 无 法 保证 其 与 
字典 原子 构成 的 子 空间 上 每 次 迭代 


的 信号 可 能 是 


就 是 说 ， 
针对 这 些 问 题 ，PATI 提 
选择 的 新 的 原子 添加 到 已 选择 


MP 算法 不 a ein 
出 了 OMP 算法 [ 


更 新 后 的 误差 向 量 与 当前 的 原子 正 交 ， 即 


因此 在 已 选 
最 优 的 。 也 


1， 该 算法 将 每 次 
原子 的 活跃 集中 ， 当 达到 要 求 的 


稀疏 度 时 停止 迭代 。 然 后 根据 最 大 
个 原子 ， 如 式 (13) 所 示 。 


PE ,3 


程度 减 小 目标 函数 值 选取 下 
1:= arg min 
| XxX Dunp lB (13) 
OMP 算法 需要 解决 的 问题 


其 中 : A 为 已 选择 的 原子 活跃 基 。 


是 | A* | 个 最 小 二 乘 问 题 , 更 新 后 的 活跃 基 、 误差 


如 式 (14) 所 示 。 


项 、 重 构 系 数 
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A:=A UI{l} 
:=(D.D,) Dix 
r:=x—D (DiD,) DX 
3.5.2 廿 优化 求解 算法 
对 于 无 约束 的 凸 优 化 问题 ， 可 以 通过 线性 规划 、 梯 度 下 降 
等 方法 凸 优化 求解 。 当 目标 函数 可 微 且 具 有 Lipschtiz07 连 续 梯 
度 函数 了 和 一 个 不 可 微 函数 462 组 成 时 ， 特 别 适 合 采用 加 速 近 


44) 


在 每 次 迭代 过 程 当中 ， 当 在 当前 点 处 线性 化 ， 求 解 转换 


mi f(a )+ Vf (a a -a)+A0() + la—a’ lg (5) 
其 中 二 次 项 为 近邻 项 ， 它 可 以 保持 更 新 后 的 了 与 线性 近似 偏差 


控制 在 一 定 范围 内 ; L 为 Lipschtiz 常数 。 式 (15) 可 以 重 写 为 
式 (16) 。 


1 ， 1 1 
min=s la—(a En 此 0 (16) 
APG 算法 求解 稀 疏 表示 问题 的 步骤 如 下 : 


关 * 1 1 
a 计算: V =4 一 一 Vf(a)，; 
L 


L 
u’ =argmin {— |u—v’ | +AQ0v’ 
b) 求 解 ; 8 ul{ 2 | |l; (7 )} ， 
SC) dv 


1\2 
9 更 新 学 习 率 ，ow 14AVL+4CD ， 


) 更 新 迭代 结果 : a”! :=w 十 人 


“一 1 
1 1-1 
Bu —u ) ; 


e)t:=ft 十 ] 若 达 到 迭代 条 件 则 退出 。 
Beck 和 Teboulle 证 明 APG 算法 的 收敛 速度 为 DG) ,证 
明 APG 算法 是 求解 稀疏 问题 的 有 效 方法 [1 引 。 


4 ”情感 识别 实验 结果 分 析 


本 文 基于 DEAP 情感 数据 库 进行 情感 分 类 。 总 共 10 位 被 
试 , 每 位 被 试 的 实验 数据 被 处 理 为 40x4x7680， 其 中 40 表示 被 
试 观看 了 40 首 音乐 视频 , 4 表示 前 额 与 额 区 的 四 个 通道 , 7 680 
表示 一 分 钟 的 采样 数据 ,根据 每 个 被 试 的 自我 评估 系统 SAM 表 
格 打 标 签 ， 可 将 情感 状态 分 为 四 个 类 别 ， 分 别 为 : 高效 价 高 唤 
醒 度 ， 高 效 价 低 唤 醒 度 ， 低 效 价 高 唤醒 度 ， 低 效 价 低 唤 醒 度 。 
基于 以 上 标签 对 情感 数据 分 别 进行 二 分 类 与 四 分 类 实验 ， 其 中 


可 以 看 出 消极 情绪 下 theta 节律 、 


9 五 种 波形 在 不 同 


青 绪 下 的 比较 


alpha 节律 、beta 节律 功 


率 谱 值 较 高 ， 而 转换 为 积极 情绪 时 delta 节律 和 gamma 节律 显 


著 增 强 。 


情感 脑 


SVM 算法 、 


二 分 类 和 四 分 类 实验 ， 


于 脑 电 信号 是 
以 其 对 非 线性 数 


居 的 处 到 


三 种 分 类 方法 均 


用 


型 的 非 线性 信和 号， 支持 向 量 机 (SVM) 
优势 以 及 较 好 的 泛 化 能 力 成 为 了 当前 
分 类 的 主流 算法 。 本 文 提 出 使 用 的 稀 玻 表示 分 类 算法 
将 与 经 典 的 SVM 算法 做 效果 比较 。 将 特征 数 
的 SPAMS 稀 玻 建 模 工 具 箱 以 及 Libsvm 工具 箱 中 ， 分 别 采 用 
通过 OMP 和 APG 算法 解码 的 稀 玻 分 类 算法 进行 
其 中 SVM 算法 选择 径 向 基 核 函数 09)， 


导入 Matlab 中 


10 折 交 叉 验 证 法 进行 验证 ， 图 10 为 三 种 分 
类 方法 的 经 过 交叉 验证 后 的 准确 率 对 比 : 


20%| 


Subpect 


三 分 类 


通过 图 10 可 知 ， 稀 玻 分 类 上 


] 分 类 


图 10 ” 稀 朴 表示 分 类 与 SVM 算法 分 类 效果 比较 


的 加 速 近邻 梯度 法 (APG) 对 


于 样本 数据 的 分 类 效果 较 好 , 二 分 类 中 最 高 准确 率 可 达到 85%， 
平均 准确 率 为 70% 左 右 ; 四 分 类 


平均 ; 


确 率 为 57.5% 左 右 。 
持 向 量 机 (SVM) 的 平均 分 类 准 


电 取 局 


确 率 可 达到 67.25%， 


而 正 交 匹配 追踪 算法 COMP ) 和 支 


确 率 要 差 些 ， 一 分 类 准 


确 率 


分 别 为 64.2% 和 64.5%, 而 四 分 类 准确 率 分 别 为 48.55% 和 44%。 


本 文 提取 了 频 域 上 的 五 个 频带 的 功率 谱 作 为 特征 , 医 


11 和 


的 分 类 准确 率 ， 


二 分 类 是 主要 效 价 维度 来 识别 被 试 的 积极 /消极 情绪 状态 。 本 文 
根据 小 波 变 换 提取 被 试 的 alpha、beta、theta、gamma、delta 节 
律 ， 并 计算 出 每 个 波形 的 功率 谱 ， 取 其 均值 作为 特征 ， 那 么 经 
过 特征 提取 后 的 数据 为 40x4x5, 其 中 5 表示 每 个 通道 中 5 种 节 
律 的 平均 功率 谱 值 。 不 同 的 波形 在 积极 情绪 与 消极 情绪 下 的 功 
率 谱 有 着 明显 的 差异 ， 如 图 9 所 示 。 


最 好 也 较为 稳定 ,APG 算法 


图 12 分 别 呈 现 了 单个 频带 作为 特 


征 时 APG 算法 和 OMP 算法 


中 ，gamma 节律 


住 确 率 能 够 达到 63.25%， 


相 较 于 其 他 四 个 波段 的 效果 
以 gamma 节律 功率 谱 作为 特征 
二 分 类 平均 准确 率 能 够 达到 64.2%， 
而 OMP 算法 在 二 分 类 时 平均 ; 


时 ， 
] 分 类 平均 准确 率 为 47%， 
] 分 类 


平均 准确 率 为 43.3%。 这 个 结果 证 实 了 gamma 节律 在 情感 识别 


中 具有 重要 的 作用 


P0] 。 
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The accuacy of APG 


The acaracy ol APG 
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二 分 类 四 分 类 
图 11 五 种 波形 分 别 在 APG 


The accuracy of OMP 


The accwacy od OMP 100% 
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二 分 类 四 分 类 
图 12 五 种 波形 分 别 在 APG 算法 下 的 准确 率 
经 过 分 析 比 较 单 个 节律 功率 谱 与 五 个 节律 功率 谱 融 合 起 来 


作为 特征 时 的 效果 ， 可 知 以 五 个 节律 功率 谱 作 为 特征 的 分 类 效 
果 要 比 单 个 节律 功率 谱 作为 特征 效果 要 好 。 


5 ”结束 语 


本 文 详细 介绍 了 基于 音乐 视频 诱导 的 DEAP 情感 数据 集 ， 
预 处 理 之 后 通过 小 波 变换 和 快速 傅 里 叶 变换 提取 其 频 域 特征 
然后 运用 通过 OMP 和 APG 算法 求解 稀 玻 编码 的 稀 玻 表示 分 类 
算法 进行 分 类 ， 并 与 SVM 算法 做 效果 比较 。 脑 电信 号 通过 稀 
疏 向 量 重 构 之 后 与 测试 样本 作 比 较 ， 通 过 宛 余 误差 进行 分 类 。 
实验 结果 表明 ， 加 速 近邻 梯度 算法 (APG) 通过 Z7 范 数 正则 近 
以 求解 以 其 快速 的 收敛 速度 在 情感 数据 集 上 有 着 较 好 的 分 类 表 
现 ， 而 贪 焚 算 法 中 的 正 交 匹配 追踪 算法 (OMP) 与 SVM 分 类 
的 效果 相差 无 几 ， 实 现 了 情感 脑 电 信号 的 分 类 。 情 感 脑 电 信号 
的 分 类 与 识别 是 一 个 值得 努力 研究 的 领域 ， 而 基于 稀 跑 表示 重 
构 信号 进行 分 类 是 对 其 有 效 的 方法 ， 值 得 再 深入 研究 。 
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